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La presente invention se rapporte a un procede de detection 
cyclique en diversite de polarisation de signaux radioelectriques numeriques 
5 cyclostationnaires concemant le domaine des radiocommunications. 

L'invention conceme en particulier la detection par deux antennes 
de diagramme de rayonnement sensible a des polarisations orthogonales de 
signaux radioelectriques. 

Dans le domaine des radiocommunications, le nombre de signaux 
10 numeriques est important et croit rapidement. Lors d'une communication 
radio un emetteur emet des signaux radioelectriques dans un ou plusieurs 
canaux frequentiels chacun caracterise par une bande frequentielle B max et 
une frequence centrale. Les signaux transmis dans ces canaux frequentiels 
sont caracterises par un certain nombre de parametres comme leur debit 

15 binaire, leur indice de modulation ou leur frequence porteuse. L'annexe A 
est un glossaire de termes utilises dans le texte. II est connu des techniques 
pour calculer la correlation spectrale ou la correlation cyclique de signaux 
cyclostationnaires pour determiner ces parametres et detecter ces signaux. 
L'annexe B regroupe un certain nombre de definitions et de proprietes 

20 concemant les signaux cyclostationnaires, les statistiques cycliques d'ordre 
2 et la correlation spectrale de tels signaux. Les systemes connus mettant 
en ceuvre ces techniques realisent cette detection a partir d'un signal recu 
sur une seule antenne. En particulier, ii est connu I'algorithme FAM, de 
("expression anglo-saxonne FFT Accumulation Method, qui a fait I'objet de 

25 publications, notamment lors du colloque intitule « 4th ASSP Workshop on 
spectrum modeling, Aug. 1988 » et en particulier sous le titre « Digital 
Implementation of spectral correlation analysers » de W.A. BROWN, H.H. 

L00MIS et sous le titre « Computational efficient algorithms for cyclic 
spectral analysers* de R.S. ROBERTS, W.A. BROWN, H.H. LOOMIS. 
30 L'algorithme FAM permet de calculer rapidement des estimateurs de 
correlation spectrale par I'intermediaire de FFT dans tout I'espace frequence 
cyclique - frequence harmonique (a,f). 

L'estimateur du premier moment E[x(f k )x(f m )*] a pour 

expression : 
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y x ,(a 0 .fo) = i;x(f k ,t)x(f m( t)*exp{-j2 J c6at} (1) 
t=l 

avec f 0 = ^ k eta 0 = f k -f m +5ct 

L'estimateur du second moment E[x(f k )x(-f m )] a pour 
expression : 

K 

5 Yx2(<*o.fo) = £x(f k ,t)x(-W)exp{-j2K8at} (2) 

t=i 

avecf 0 = ^-y^eta 0 =f k +f m +5a 

Les signaux x(f,t) sont calcules par une batterie de filtres de 
bande Be** et de frequence centrale notee f k , f m . Le signal x(f,t) a done une 
bande valant Bear*. Pour se rapprocher de la representation frequentielle 

10 exacte du signal x(t), il faut que soit tres inferieur a la bande B du 
signal x(t). Comme x(f) est calcule dans une bande non nulle ce signal x(f) a 
une evolution au cours du temps et depend alors de f et t : x(f,t). Dans une 
premiere etape, les signaux x(f k) t) et x(f m ,t)* sont intercorreles apres avoir ete 
mis en bande de base. Ces signaux etant dans une bande B„n«i non nulle, le 

15 signal obtenu z(t) evolue au cours du temps dans une bande de largeur 
2xB«rm. Etant donne que z(t) est dans une bande non nulle et que le pic de 
correlation en frequence cyclique n'est pas forcement positionne en f^fm, il 
est necessaire de realiser une transformee de Fourier de z(t) pour obtenir le 
decalage 5a de ce pic de correlation. 

20 Pour permettre un calcul rapide des estimateurs de correlation 

spectrale, et y^, I'algorithme FAM se deroule en plusieurs etapes. 

Dans une premiere etape, les signaux x(fic,t) et x(-f m ,t) sont 
calcules par un systeme de FFT glissantes. L'ecart Aa = f k+1 -f k entre les 

frequences et f k est constant et la frequence d'echantillonnage des 
Fe 

25 signaux x(f,t) vaut Fe^tf = — , avec RE le recouvrement des fenetres 

glissantes en nombre d'echantillons et Fe la frequence d'echantillonnage de 
x(t). 

Dans une deuxieme etape, est effectue le calcul du signal 
z(t) = x(f k ,t)xx(f m ,t)* pour le premier moment et du signal 
30 z(t) = x(f k , t) x x(-f m , t) pour le second moment. 
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Dans une troisieme etape, est realisee une FFT du signal z(t) sur 
un certain nombre d'echantillons. 

Sur le premier moment pour calculer la correlation spectrale en 
f _ ( f k + f m) et pQur des a _ oompris entre a min =f k -f m -A^- et 

5 a max = f k ~ f m s " r K echantillons, il faut realiser une FFT du signal 
z(t) = x(f kl t) x x(f m ,t)* sur K echantillons. Les K echantillons frequentiels du 
signal z(a«) resultant representent la correlation spectrale en f = ^ k * fm ) 
avec des a* compris entre ct min =f k -f m - Fe< ^ nal et 
a max = f k ~ f m + Fe ^ nal . Comme Fe^ est plus grand que Act seules sont 
10 retenues les frequences cycliques a* contenues dans I'intervalle 
(fj< -f m - A^j —(f K -f m + A~j . Toutes les sorties x(f k ,t) et x(f m ,t) du banc 

de filtres sont corr&ees pour obtenir les points de la correlation spectrale 
dans la zone d'interet du premier moment c'est-a-dire pour a ( B max . 

Sur le second moment, le calcul de la correlation spectrale est 
1 5 limite en frequence harmonique par B^. 

L'algorithme FAM expose prec6demment ne s'applique qu'a des 
signaux issus d'une seule antenne. Or, dans le domaine des 
radiocommunications les signaux se propagent avec une certaine 
polarisation et toute antenne a un gain qui depend de la polarisation des 
signaux recus. II en requite qu'une seule antenne filtre « en polarisation » 
les sources et peut meme completement annuler les signaux recus ce qui 
interdit toute possibility de detecter les sources. Et ceci quel que soit 
l'algorithme utilise, en particulier en utilisant un algorithme FAM. 

Le but de I'invention est de proposer une solution a ce probleme. 

A cet effet, I'invention a pour objet un procede de detection 
cyclique en diversite de polarisation de signaux radioelectriques numeriques 
cyclostationnaires de frequence d'echantillonnage Fe transmis dans un 
canal frequentiel de bande B,,,,, et recus sur un reseau de N superieur ou 
egal a deux antennes dont leur diagramme de rayonnement presente un 
maximum maximorum de sensibilite caracterise en ce qu'il consiste, pour 
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tout couple d'antennes, a acquerir sur une duree d'observation T et dans 
une bande d'acquisition B««, les signaux numeriques en sortie des antennes, 
a calculer, pour chaque frequence cyclique d'un espace frequence 
cyclique - frequence harmonique determine limite par la bande du canal 
5 frequentiel et la bande d'acquisition B,*, des signaux numeriques acquis, 
une matrice de covariance cyclique sur un premier moment, respectivement 
sur un second moment, d'une correlation spectrale des signaux numeriques 
acquis pour detecter des pics de correlation spectrale en comparant un 
rapport de vraisemblance determine a partir de la matrice de covariance 
1 0 cyclique avec un seuil de detection determine d'une maniere statistique. 

L'invention consiste en un test de detection cyclique de signaux 
radioelectriques recus par tout couple d'antennes d'un reseau de 
N antennes dont leur diagramme de rayonnement presente un maximum 
maximorum de sensibilite. 
15 Dans un premier mode de mise en ceuvre du precede, le reseau 

comprend N=2 antennes dont les maximums maximorum de sensibilite 
pointent dans des directions orthogonales. 

Dans un deuxieme mode de mise en ceuvre du procede, le reseau 
comprend N=2 antennes dont les maximums maximorum de sensibilite 
20 pointent dans une meme direction. 

L'invention a pour principal avantage de realiser un test de 
detection rapide de signaux radioelectriques presentant une polarisation 
quelconque en calculant un rapport de vraisemblance suivant une relation 
simple. Dans le cas d'un reseau a deux antennes et en absence de source, 
25 le rapport de vraisemblance suit une loi statistique du chi_2 a huit degres de 
liberte. La connaissance de la loi de probability du rapport de vraisemblance 
en I'absence de source permet de calculer un seuil de detection avec une 
certaine probability de fausse alarme pour juger de la presence d'une 
source en un endroit determine de I'espace frequence cyclique - frequence 
30 harmonique. 

Le test de detection cyclique du rapport de vraisemblance est 
d'une part sur deux antennes et d'autre part statistique. En consequence, la 
detection ne se fait pas avec un seuil de detection empirique. 

Le procede selon l'invention a aussi pour avantage de penmettre 
35 la differentiation de signaux radioelectriques presents dans un meme canal 
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frequentiel mais transmis avec des debits differents. Le procede selon 
i'invention permet avantageusement de separer les signaux radioelectriques 
selon leur frequence de modulation mais aussi selon leur debit de 
transmission. 

5 D'autres avantages et caracteristiques de I'invention apparaitront 

dans la description qui suit faite en regard des figures annexees qui 
represented : 

-la figure 1, le procede selon I'invention sous la forme d'un 
organigramme, 

10 - la figure 2, une projection sur I'axe a des frequences cycliques 

du rapport de vraisemblance dans le cas de trois sources correspondant a 
des canaux frequentiels differents et transmises a un meme debit de 
25 kbits/s, 

- les figures 3 et 4, le resultat du test de detection dans I'espace 
15 frequence cyclique - frequence harmonique effectue sur le rapport de 

vraisemblance represents sur la figure 2, 

-la figure 5, un trace dans I'espace frequence cyclique - 
frequence harmonique delimite par la figure 4 du niveau de la correlation 
spectrale du signal aux points (a*,fk) ou a ete detectee une source, 

20 - la figure 6, une projection sur I'axe a des frequences cycliques 

du rapport de vraisemblance dans le cas de deux sources correspondant a 
un meme canal frequentiel et transmises a des debits differents, 

-la figure 7, le resultat du test de detection dans I'espace 
frequence cyclique - frequence harmonique effectue sur le rapport de 

25 vraisemblance represents sur la figure 6, 

- la figure 8, un algorithme sur le premier moment d'un procede 
selon I'invention, 

- la figure 9, un algorithme sur le second moment d'un procede 
selon I'invention, 

30 - la figure 10, la correlation spectrale du premier moment d'un 

signal avec une modulation lineaire BPSK filtree Nyquist, 

- la figure 1 1 , la correlation spectrale du second moment d'un 
signal avec une modulation lineaire BPSK filtree Nyquist, et 
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- les figures 12a et 12b, les zones d'interets frequence cyclique - 
frequence harmonique des correlations spectrales sur les premier et 
second moments d'un signal de bande finie B. 

Sur les differentes figures, les elements homologues sont 
5 represents avec une meme reference. 

Le procede selon I'invention se deroule de la facon representee a 
la figure 1 suivant quatre etapes principales. Dans une premiere etape, 
representee en 1 a la figure 1, il consiste a recevoir des signaux 
radioelectriques numeriques de frequence d'echantillonnage Fe, emis par 
10 une source dans un canal frequentiel de transmission de bande B max , par 
un reseau de N antennes A1, ... AN avec N superieur ou egal a 2, N=2 dans 
I'exemple represents sur la figure 1, et a acquerir sur une duree 
d'observation T et dans une bande d'acquisition B acq les signaux de sortie 

des antennes pour les ranger dans des vecteurs x(f k< t). L'acquisition 1 des 
15 signaux numeriques en sortie des antennes A1, AN consiste a filtrer chaque 
signal numerique recu sur une antenne par paquets de K echantillons par 
une batterie de M filtres de bande elementaire B^nai, dont les sorties sont 
sous-echantillonnees d'un facteur RE, et de frequence centrale donnee 
notee f ta f,„ avec k et m variant entre 1 et M, K est donne par la relation : 
TxFe 

20 K = . La batterie de M filtres est telle que I'ensemble des frequences 

RE 

centrales realise un decoupage regulier d'une bande d'acquisition B acq de 
maniere a verifier la relation B acq ^MxBcanai. Le filtrage est realise par la 
mise en ceuvre de FFT suivant des techniques connues de I'homme du 
metier. Apres filtrage chaque signal numerique filtre est range dans des 
25 vecteurs xft.*) dont les composantes x^f^.t) x 2 (fk,t) .... x N (f k ,t) 

correspondent aux signaux des differentes antennes A1, AN filtres par un 
meme filtre de frequence centrale determinee f*. Les vecteurs x(f k ,t) 
permettent lors de I'etape 2 de calculer des matrices de covariance cyclique 
Utah) sur un premier moment, respectivement sur un second moment, d'une 
30 correlation spectrale des signaux numeriques recus. Lors d'une troisieme 
etape, representee en 3, un rapport de vraisemblance V(a*) est determine a 
partir des matrices de covariance U(at,) et est compare lors de I'etape 4 a un 
seuil determine d'une maniere statistique par une loi du chi_2. La 
comparaison permet de determiner les doublets frequence 
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harmonique - frequence cyclique (f 0l oto) pour lesquels une source au moins 
a ete detectee. 

Lors d'une communication radio, en provenance d'un satellite par 
exemple, ia polarisation des signaux radioelectriques evolue en fonction du 
5 milieu de propagation. Le reseau d'antennes utilise comprend un nombre N 
superieur ou egal a deux antennes. Ces antennes ont pour propriete de 
presenter un diagramme de rayonnement presentant un maximum 
maximorum de sensibilite dans une direction donnee. 

Dans un premier mode de mise en ceuvre du precede, le reseau 
10 comprend un nombre N d'antennes egal a 2 dont les maximums maximorum 
de sensibilite porntent dans des directions orthogonales. Dans d'autres 
modes de realisation de I'invention, le procede peut etre mis en oeuvre pour 
tout couple d'antennes d'un reseau de N antennes ayant un diagramme de 
rayonnement presentant un maximum maximorum de sensibilite dans une 
15 direction donnee. 

En presence d'une source, correspondent par exemple a 
remission d'un satellite, de direction d'arrivee azimut-site (6o,Ao) et de 
polarisation P 0< le signal recu sur I'antenne « n » est le suivant : 

x n (t) = s(t)a n (e 0 ,A 0l Po)+b n (t) (3) 
20 ou a n (9o,Ao,Po) est la reponse de I'antenne « n » a la direction 

(6o,Ao) et a la polarisation P 0 , 

b„(t) le bruit additifet 

s(t) le signal temporel de la modulation du signal emis par la 

source. 

25 Le signal recu sur I'ensemble des N antennes est un signal 

vecteur de dimension N dont I'expression est la suivante : 
x(t) = s(t)a(eo.Ao,Po) + b(t) (4) 

a(0 o , A 0 .P 0 ) = [ai(e 0 . A 0 .P 0 ) - aN(eo. Ao.Po)]* . 
30 b(t) = [b 1 (t)-b N (t)] t . 

Le vecteur a(6o,Ao,Po) est appele vecteur directeur. La 
polarisation P 0 se decompose en une composante horizontale P H0 et une 
composante verticale Pvo- Le vecteur directeur a la caracteristique de se 
decomposer lineairement suivant ces deux composantes, soit : 
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15 



20 



a(e 0 ,Ao,Po) = Pho aH(e 0 ,Ao)+Pvo M Q o,*o) 

a H (e 0 ,A 0 ) : reponse du reseau d'antenne en polarisation 

horizontale pour une direction (G 0 ,Ao), 

a v (e 0 ,A 0 ) : reponse du reseau d'antenne en polarisation 

verticale pour une direction (9 0 ,Aa). 

Apres detection, une goniom6trie est generalement mise en 
ceuvre pour localiser les sources detectees. Etant donnee la decomposition 
de la polarisation P 0 , la goniometrie est realisee a partir de deux tables de 
calibration : la premiere pour une polarisation horizontale et la seconde pour 
une polarisation verticale. 

Sur le premier et le second moment de la correlation spectrale les 
matrices de covariance cyclique sont definies suivant les relations : 



pour le 1* moment : Rx(a,f) = E 



pour le 2 nd moment : C^aJ) = E 



(6) 



Les matrices de covariance cyclique font intervenir toutes les 
intercorrelations prises 2 a 2 entre les antennes. Ainsi, le terme de la ni*"" 
ligne et de la na 1 *"" colonne des matrices de covariance cyclique represente 
I'intercorrelation cyclique entre les antennes ni et n 2 . 

pour le 1 W moment : 



jr le 2 nd moment : 

^wW f+ fMf-')] 



pour le 2' 
C 



(7) 



(8) 



En supposant que le bruit soit stationnaire, les matrices de 
covariance cyclique s'expriment de la facon suivante en fonction du vecteur 
25 directeur a(9o, Ao.Po) et de la correlation spectrale du signal s(t) : 

pour le 1" moment : 

Rx(a,f) = Ysl(a,f) a(8 0l Ao.PoJafeo.Ao.Por (9) 
pour le 2 nd moment : 

C x (a,f) = y S 2(a,f) ^eo.Ao.PoJaCBo.Ao.Po) 1 (10) 
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Ainsi, lorsque le signal s(t) a une correlation spectrale sur le 
premier ou le second moment non nulle, la matrice de covariance cyclique 
associee est de rang 1. En consequence, les matrices de covariance 
cyclique ont les memes proprietes cycliques que la correlation spectrale du 
5 signal a goniometrer. 

Dans la realite ces matrices de covariance cyclique sont estimees 
sur BT echantillons independants de la meme facon que les correlations 
spectrales et valent : 

pour le 1* moment: 

10 R*(fk -fm ^^^'^fkM^T^-r^^) ( 11 > 

pour le 2 nd moment : 

C x (fk +*n +&t.fi L ~i) = S^.W-fm-t) 1 exp{-j27c6a t} (12) 

Ainsi lors du traitement, les triplets (fk,f m ,5a) detectes 
cycliquement sont stockes puis reutilises pour integrer les matrices cycliques 
15 et realiser la goniometrie. En consequence, il faut stocker les signaux 
vectoriels x(f p ,t) et x(-f p ,t) dans un fichier. II faut remarquer que les 
echantillons x^.t) et x(fk,t+l) ne sont pas forcement independants, c'est 

pourquoi K est superieur ou egal a BT. La relation entre K et BT depend de 

Fe 

la frequence d'echanti I tonnage Fe^, Fe^a) = — , et de la bande Bca* 

RE 

20 des signaux x(f k ,t) comme le montre ('expression suivante : 
BT= Bca na ,xK (13) 
F e canal 

Pour connaTtre les couples (<Xn,f„) pour lesquels est detectee au 
minimum une source, un test de detection cyclique selon I'invention sur deux 
antennes du nombre de sources multicapteurs est mis en ceuvre. Ce test de 
25 detection cyclique depend du produit BT, ou B = B cana) correspond a la 
1^ 

bande des signaux x(f.t) et T = correspond au temps d'integration 

^ e canal 

des matrices de covariance qui coincide avec la duree d'observation. 

Ce test de detection cyclique a pour objectif de compter le nombre 
de sources presentes en meme temps pour le point frequence 
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cyclique - frequence harmonique (ao,fo). En presence de M sources en 
(ao,fo), le modele du signal devient : 

M % 
x(t)= Z s m(t)a(e mi A mi P m ) + Zs m .(t)a(e m . I A m .P m .) + b{t) (14) 

m=1 m'=M+1 
D'apres ce modele M T sources sont recues dont M ont de 
5 I'energie pour leur correlation spectrale en (<xo,fo). Suivant ce modele, en 
classlque (utilisation de la matrice de covariance) M T sources doivent etre 
goniometrees alors qu'en cyclique (utilisation de la matrice de covariance 
cyclique) seulement M sources doivent etre goniometrees. En effet, les 
matrices de covariance cyclique ont la caracteristique d'isoler une partie des 
10 sources. 

Ce test de detection cyclique consiste a estimer ie nombre de 
valeurs singulieres non nulles d'une matrice d'intercorrelation entre les 
signaux & et yj permettant ainsi d'estimer le rang de cette matrice. Suivant 
que ie premier ou le second moment est utilise^ ces signaux & et y, valent : 
15 1"moment : x t = xff 0 +^-,t] et y t = xff 0 -^-.t) (15) 

2" d moment : x t = x(f 0 + ^-,t) et y t = x^-fo.t) (16) 

Cette intercorrelation se fait sur K echantillons qui correspondent 
a BT echantillons statistiquement independants : 
1 K 

R xy = T7X*y t + (17) 
K t=i 

20 Le rang de Rxy est estime par le rang de sa forme normalisee R. 

Sachant que les signaux temporels & et y, sont des processus gaussiens, en 
presence de Mo sources cycliques le rapport de vraisemblance des N-M 0 
plus faibles valeurs propres de R suit une Joi du chi_2 a 2(N-M 0 ) 2 degres 
de liberte : 



25 



Vcydk,ue([M = Mo]/M 0 )=-2BTIr( fl*m } (18) 
avec R = I-UU + 



o6U = Ri«R w R;« 
M<n : valeurs propres de R 
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N : nombre d'antennes 

Le test de detection cyclique consiste a faire un test de detection 
presque classique sur la matrice hermitienne normalisee R. Pour cette 
matrice les valeurs propres de baiit sont exactement egales a 1. C'est 
5 pourquoi dans ce rapport de vraisemblance, le logarithme de la valeur 
moyenne du bruit n'est pas retranche. La seule difference par rapport au 
classique est le nombre de degres de liberte : 2(N-M 0 ) 2 au lieu de 
(N-M 0 ) 2 -1. 

La connaissance des statistiques de V cycJ j que (M 0 /M 0 ) permet 
10 de fixer un seuil p M pour lequel la probability d'avoir strictement plus de M 
sources est proche de 1 (pd~1). Le test est le suivant : 

• si V(M/M 0 ))Pm- le nombre de sources M 0 presentes est 
superieur a M, 

• si V(M/M 0 )(Pm. Ie nombre de sources Mo presentes est 
15 inferieur ou egal a M. 

Pour determiner le nombre des sources Mo, le test teste d'abord la 
presence de M=0 source puis de M=1 source jusqu'a ce que V(M/M 0 ) soit 

inferieur au seuil 0m. Le nombre de sources vaut alors Mo=M. 

Le procede selon I'invention s'applique pour tout couple 

20 d'antennes, soit N=2. Le rapport de vraisemblance defini ci-dessus par la 
relation (18) est calcuie avec les antennes c 1 » et « 2 », dans tout I'espace 
(ao.fo) frequence cyclique - frequence harmonique puis un seuil de detection 
d'au moins une source est fixe. En I'absence de source le rapport de 
vraisemblance du test de detection cyclique suit une loi du chi_2 a 

25 2 x 2 2 = 8 degres de liberte et devient : 

V xd 2 C ydiq U e(a 0 .fo) = -2BT lnTdet[l-U xd (ao,fo) U xri (a 0 ,f 0 ) + ]] (20) 
ou sur le 1* moment : U^(a 0 ,f 0 ) = R z (a 0 ,fo) (21) 
ou sur le 2 nd moment : U x2 (a 0 ,fo) = C z (a 0 ,fo) (22) 
(d prenant la valeur 1 pour le premier moment, et la valeur 2 pour 
30 le second moment) 

ou z{f,t) = R x (f)- 1/2 x(f,t) (23) 
ou detrj designe le determinant 
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ou R x (fo) la matrice de covariance a la frequence f 0 du signal 
x(t) = [x 1 (t)x 2 (t)] t 

et R z (a 0 ,fo) et C z (a5,fo) las matrices de covariance cyclique 
du signal z(t) = [z 1 (t)2 2 (t)] t . 

La connaissance de la loi de probabilite du rapport de 
vraisemblance en I'absence de source permet de calculer le seuil de 
detection avec une certaine probabilite de fausse alarme : 

Lorsque N=2, la matrice R x (f)~ V2 a par exemple I'expression 



litterale suivante : 
10 R x (f)" 



)/ R *( f )(l.l) 




(24) 



Cette derniere expression permet de calculer facilement le rapport 
de vraisemblance V xd 2 C yciique(a,f) dans tout I'espace frequence cyclique - 
frequence harmonique. Le calcul de V xd 2 C yclique(ao.fo) est simplifie par la 
connaissance de ('expression litterale du determinant de la matrice 2x2 : 
15 I - U xt j(a,f) U xc j(a,f) + . Le test de detection cyclique est le suivant : 

Si V xd 2 cyC |jque(a 0l f 0 )>seun[x 2 (8>p fa ] (25) 

alors il y a une source en (oco,fo). 

Le seuil de detection est dans ce cas determine a partir de la loi 
du chi_2 avec 8 degres de liberie pour une certaine probabilite de fausse 
20 alarme p* la plus faible possible. Le test de detection cyclique selon 
invention est illustre par des exemples issus de simulation et representes 
par les figures 2 a 7. 

Les figures 2, 3, 4 at 5 illustrent le cas de sources de meme debit, 
25 kbits/s mais dont les frequences de modulation sont differentes. 
25 Les figures 6 et 7 illustrent le cas de sources de debits differents 

mais qui se trouvent dans un meme canal frequentiel. 

Sur la figure 2 est represents le rapport de vraisemblance, sur le 
premier moment, en presence de trois MSK de meme debit egal a 25 kbits/s. 
Le temps d'integration est de BT=150 echantillons independants. 
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Pour huit degres de liberte !e seuil de detection est situe en -7. 
Tous les pics de correlation dont le niveau est superieur a ce seuil sont 
detectes. Le test de detection cyclique donne un resultat illustre sur les 
figures 3 et 4, dans I'espace frequence cyclique - frequence harmonique. Le 
5 test de detection cyclique permet de detecter les trois sources MSK 

Un pic de correlation est obtenu tous les 25 kHz car les 3 MSK 
ont un debit de 25 kbits/s. Pour confirmer qu'il s'agit bien de MSK, un trace 
du niveau de correlation spectrale du signal aux points (ct*,f k ) ou a ete 
detectee une source peut etre realise. II est illustre par la figure 5. 
10 La figure 6 represente le rapport de vraisemblance sur I'axe des 

frequences cycliques a sur le premier moment. Pour huit degres de liberte, 
le seuil de detection est situe en -7. Tous les pics de correlation dont le 
niveau est superieur a ce seuil sont detectes. Ceci permet de detecter un pic 
de correlation a a = 16 kHz et un pic de correlation a a = 25 kHz. Le test de 
15 detection cyclique donne un resultat illustre sur la figure 7 dans I'espace 
frequence cyclique - frequence harmonique. Le test de detection cyclique 
permet de detecter les deux sources situees dans le meme canal frequentiel 
et de les differencier par leur debit. 

Finalement, le procede selon I'invention peut etre schematise 
20 d'une maniere condensee illustree par la figure 8 pour le premier moment et 
par la figure 9 pour le second moment. Les vecteurs signaux x(f k ,t) et 
xftn.t) . obtenus lors de I'acquisition par tout couple de N = 2 antennes, sont 
blanchis 5 avec les matrices R x (f k )" 1/2 et R x (f m ) _1/2 donnees par la 
relation (24) pour obtenir les vecteurs signaux z(f k ,t) et z(f m ,t). Ces 
25 vecteurs signaux permettent de calculer les K matrices U(t) = z(f k , t) . z(f m , t) + 
de dimension 2x2 dont les composantes sont notees u^. 

Une FFT 6 est realisee pour chaque composante u y sur K points 
en prenant les K composantes u (j des K matrices U(t). Etant donne que les 
matrices U(t) sont de dimension 2x2, ceci conduit a realiser quatre FFT sur 
30 K points. 

Les K matrices UK) resultantes sont de meme dimension 2x2 
que les matrices U(t). 

Les K matrices U(cck) permettent le calcul 3 du rapport de 
vraisemblance V(a*) suivant la relation : 
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V(ct k ) = -2 . BT In [det[l - U(a k ) U(a k ) + J 

Cette expression simple permet de calculer rapidement le rapport 
de vraisemblance a la frequence harmonique f = ^ * fm ^ et aux frequences 

cycfiques a* comprises entre a,,^ = f k - f m - et a max = f k - f m + a|- . 

5 Les doublets de frequences (f kl f m ) sont tels qu'ils permettent de 

calculer le rapport de vraisemblance pour les frequences cycliques 
inferieures a ow 

Le test de detection cyclique, suivant la relation (25), est realise 
sur chacun des points du rapport de vraisemblance calculus precedemment. 

10 Ceci permei de deduire les doublets a k j pour lesquels au moins 

une source a ete detectee. 

Sur le second moment, figure 9, le principe de I'algorithme est le 
meme que celui pour le premier moment en correlant les vecteurs z(f k .t) et 
z(-f m ,t). Le calcul 3 du rapport de vraisemblance est dans ce cas limite sur 
1 5 les frequences harmoniques par Bm«. 
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Annexe A 
GLOSSAIRE 

5 A.1. Conventions standards : 

x* : conjuguee de x 

x l : transposee du vecteur x 

x* : transposee conjuguee du vecteur x 

X* : transposee de la matrice X 
10 X* : transposee conjuguee de la matrice X 

X' 1 : inverse de X 

E{.] : esperance mathematique 

x, X : valeurs estimees 

x* y : convolution des signaux x et y 
15 ® : produit de Kronecker 

BPSK : bi-phase shift keying : modulation lineaire a deux etats de 

phase 

MSK : minimum shift keying : modulation par deplacement de 
frequence a indice de modulation 0,5 et a phase continue 
20 FSK : frequency shift keying : modulation par deplacement de 

frequence 

FFT : Fast Fourier Transform : transformee de Fourier rapide 
In : logarithme neperien 

25 A.2. Constantes entieres : 

N :nombre d'antennes 

M T : nombre de sources a goniometrer suivant les methodes 

classiques 

M : nombre de sources a goniometrer suivant I'invention 
30 B : bande frequentielle 

Bcami : bande frequentielle elementaire d'un fiitre d'une batterie de 

filtres. 

Bmax : bande frequentielle d'un canal frequentiel do transmission. 

35 A.3. Variables scalaires : 

x„(.) :observation sur I'antenne « n » 
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s() : signal temporal de la modulation du signal emis par la source 
b„() : bruit additif sur I'antenne « n » 
Fe : frequence d'echantillonnage de x(t) 

a„(8o,Ao,Po) : reponse de I'antenne « n » a la direction (9 0l Ao) et a 
5 la polarisation P 0 

6 0 : azimut de la direction d'arrivee du signal emis par la source 
Ao : site de la direction d'arrivee du signal emis par la source 
S m () : signal temporal de la modulation du signal emis par la 
source « m » 

1 0 f p , f 0l f k , f m , f : frequences harmoniques 

a, ao : frequences cycliques 

RE : recouvrement des fenetres glissantes en nombre 
d'echantillons 

M«, : valeur propre de R 

15 

A.4. Variables vectorielles et matricielles : 

x : signal vecteur de dimension N correspondant au signal recu 
sur les N antennes 

I : matrice identite 
20 R x : matrice de covariance du vecteur x 

R y : matrice de covariance du vecteur y. 

a(8 0 ,Ao,Po) : vecteur directeur de dimension N correspondant a la 
reponse des N antennes a la direction (9 0 ,Ao) et a la polarisation P 0 
Po : polarisation du signal s(t) 
25 R^ : matrice d'intercorralation entre les vecteurs x, et & 

R : matrice hermitienne normalisee de la matrice Rxy 



2774217 



17 



Annexe B 

Signaux cyclostationnaires recus sur un capteur. 
Statistiques cycliques d'ordre 2 et correlation spectrale. 

5 

Un signal est non stationnaire lorsque sa distribution statistique 
depend du temps. Parmi cette classe, les signaux cyclostationnaires ont des 
statistiques evoluant periodiquement au cours du temps. Dans le cas de 
signaux numeriques (BPSK, FSK, etc.) ces periodicites dependent en 
10 particulier du debit de la modulation et de la frequence porteuse. Les 
statistiques d'ordre 2 d'un signal x(t) sont entierement caracterisees par le 
premier moment r^t.t) et le second moment c^t.t) tel que : 

1* moment: r x (tt) = E< 



20 



2 nd moment: c x (t,t) = E 

15 Un signal est stationnaire a I'ordre 2 lorsque les moments 

associes sont independants de t, soit : 
r x (t.x) = r x o(x) 
etc x (tT) = c x0 (t) 

Un signal est cyclostationnaire a I'ordre 2 lorsque ces moments 
sont periodiques en t, soit : 
r x (t.T) = r x (t + T 1 ,T) 

etc x (t,T) = c x (t+T 2 ,T) 

II faut remarquer que les periodicites Ti et T 2 des premier et 
second moments ne sont pas forcemeat egales. Ces deux moments peuvent 

alors s'exprimer suivant une serie de Fourier, soit : 

r x (t, t) = ^ r x (na-j, x) exp(j27t na^) avec «i = 

c x(t x ) = Z c x(na2. exp(j2re na 2 t) avec a 2 = 

n '2 

25 Les coefficients de Fourier r K (a,x) et c„(a,T) constituent les 

fonctions d'autocorrelation cyclique. Us admettent des pics de correlation 
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15 



20 



pour les frequences cycliques a multiples de respectivement a-j = — et 

M 

a 2 =~- et sont les transformers de Fourier en « t » des moments r x (t,t) et 
T 2 

Cx(t,T), soit : 



r x (a,T)=lim 



c x (a 1 t)=lim 



1 Jr x (t,t)exp{-i27tat}dt 

T 

1 Jc x (t,T)exp{-j27tat}dt 

1 T 



Les fonctions d'autocorrelation cydique r^a,*) et Cx(a,x) admettent 
pour transforms de Fourier en x les correlations spectrales yi(a,f) et y 2 (a,f), 
tel que : 

YxlM) = f r x (a,t) exp{-j2* ft) dt 
Yx2(a,f) = Jc x (a,T)exp{-j27cft}dT 

X 

Si les fonctions r x (a,t) et Cx(a,x) admettent des pics de correlation 
pour les frequences cycliques a multiples de respectivement et a 2 , il en 
sera de meme pour y x i(a,f) et yx2(a,f). Les fonctions yxi(a,f) et yx2(a,f) sont des 
correlations spectrales car elles ont pour expression : 



1* moment: y x1 (a l f) = E 



2 nd moment: Yx2(a,f) = E 



<-f)<f-0] 



II faut remarquer que sur le premier moment en a=0, nous 
obtenons la densite spectrale du signal, c'est-a-dire son spectre de 
frequence. 

Les signaux non stationnaires ont des moments r x (t,t) et Cx{t,x) 
dependant du temps « t », il en resulte que les correlations spectrales y*i(a,f) 
et yx2(a,f) sont non nulles en dehors de a=0 sur le premier moment. 
Contrairement aux signaux cyclostationnaires il n'y a pas de pic de 
correlation mais des continuums d'energie sur I'axe frequence cyclique a. 
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Dans le cas d'un signal de bande finie B, I'espace de calcul du 
premier moment et du deuxieme moment petit etre limite a un espace 
determine. Soit I'exemple d'une modulation lineaire et en particulier d'une 
modulation BPSK filtree Nyquist. Cette BPSK a un debit normalise de 
00=0,25 avec une frequence porteuse normalisee de fo=0,1. 

La figure 10 represente la correlation spectrale du premier 



moment : 



1* moment (BPSK) : Yx1 (a,f) = E 



f + 7 f" 



10 



La figure 10 illustre le fait qu'il y a presence d'energie pour les 
frequences cyciiques a multiples du debit eta D'autre part, le spectre en 
frequence harmonique est centre en f=f 0 . En consequence, lors de la 
detection des couples (o^fO presentant de I'energie, il faut scruter les 
frequences harmoniques presentes dans la bande d'acquisition du signal. II 
faut de plus remarquer que la correlation spectrale est symelrique autour de 
15 I 'axe <x=0. 

La figure 11 represente la correlation spectrale du second 

2 nd moment : yx2 (a,f) = E^f + 1] >{|-f]j 

Avec le second moment, il y a presence d'energie pour les 
20 frequences cyciiques a egales a ocok+2f 0 (k : entier). Ainsi les pics de 
correlation en frequence cyclique vont dependre a la fois du debit et de la 
frequence porteuse. Pour la detection d'energie, il faut done scruter les 
frequences cyciiques entre -B«, et B^ (B^ est la bande d'acquisition du 
signal). D'autre part, le spectre en frequence harmonique est centre en f=0. 
25 En consequence, lors de la detection des couples(ot kl f lt ) presentant de 
I'energie, il faut scruter les frequences harmoniques presentes dans la 
bande B des signaux qui peut etre tres inferieure a la bande d'acquisition 
Bacq. II faut remarquer que la bande du signal en frequence harmonique est 
de I'ordre de grandeur du debit. 
30 Cet exemple illustre clairement que pour un signal de bande finie 

B et de frequence porteuse f, (ou f 2 ), les correlations spectrales sur les 
premier et second moments sont non nulles dans des zones frequence 
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cyclique - frequence harmonique dependantes de f^ou f 2 ) et B. Les figures 
12a et 12b represented schematiquement ces zones interessantes. 

La bande du signal emis n'est pas toujours connue de la 
reception, par contre le canal frequentiel de transmission utilise est connu 
5 par sa bande B max , caracteristique du filtre d'emission de I'emetteur, et par 
sa frequence centrale. D'apres ces remarques, de I'energie est presente de 
facon significative sur les premier et second moments de la correlation 
spectrale lorsque ; 

— B B 

10 1"moment: 0<a<B max et — |^.^f^-^. 

2-moment: -B acq <a<B acq et ^.f.?f 

En consequence, il faut realiser la detection cyclique uniquement 
dans ces zones frequence cyclique - frequence harmonique . 
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REVENDIC ATION S 

1 . Procede de detection cyclique en diversite de polarisation de 
5 signaux radioelectriques numeriques cyclostationnaires de frequence 
d'echantillonnage Fe transmis dans un canal frequentiel de bande B™x et 
recus sur un reseau de N superieur ou egal a deux antennes (A1.AN) dont 
leur diagramme de rayonnement presents un maximum maximorum de 
sensibilite caracterise en ce qu'il consiste, pour tout couple d'antennes 

10 (A1.AN), a acquerir (1) sur une duree d'observation T et dans une bande 
d'acquisition les signaux numeriques en sortie des antennes (A1.AN), a 
calculer (2), pour chaque frequence cyclique d'un espace frequence 
cyclique - frequence harmonique determine limite par la bande du canal 
frequentiel et la bande d'acquisition B^ des signaux numeriques acquis, 

15 une matrice de covariance cyclique sur un premier moment, respectivement 
sur un second moment, d'une correlation spectrale des signaux numeriques 
acquis pour detecter (4) des pics de correlation spectrale en comparant (4) 
un rapport de vraisemblance (V(<xk)) determine (3) a partir de la matrice de 
covariance cyclique (U(ak)) avec un seuil de detection (p M ) determine d'une 

20 maniere statistique. 

2, Procede selon la revendication 1, caracterise en ce que 
I'acquisition (1) des signaux numeriques consiste a filtrer, par une batterie 
de M fiitres de bande elementaire B^ai dont les sorties sont sous- 
25 echantillonnees d'un facteur RE et de frequence centrale donnee (f k ) telle 
que I'ensemble des frequences centrales (f k ) realise un decoupage regulier 
de la bande d'acquisition B^ de maniere a verifier la relation 
B a cq ^MxBcanai, chaque signal radioelectrique numerique repu sur une 

antenne par paquets de K echantillons avec K donne par la relation 

30 K = et a ranger chaque signal numerique filtre dans des vecteurs 

Rc 

(xft.t)) dont les composantes (x 1 (f k ,t) x 2 (fk,t), .... x N (f k ,t)) correspondent 
aux signaux des differentes antennes (A1.AN) fiitres par un meme filtre de 
frequence centrale determinee (f k ). 
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3. Procede selon la revendication 2, caracterise en ce que le 
calcul (2) d'une matrice de covariance cyclique (U(ak)) sur le premier 
moment consiste a choisir deux vecteurs (x(fk,t) x(f m ,t)) de frequence f k et 
f m , a les blanchir (5) pour obtenir deux nouveaux vecteurs (z (f^ . z(f m ,t)), 

5 a calculer K matrices U(t) donnees par la relation U(t) = z(f kl t) x z(f m ,t) + 
correspondant a la duree d'observation T, a calculer (6) chaque composante 
UiKoK) de la matrice de covariance cyclique (U(ak)) en realisant une 
transformee de Fourier rapide sur K points en prenant les K points de meme 
indice ij des K matrices U(t). 

10 

4. Procede selon Tune quelconque des revendications 2 et 3, 
caracterise en ce que le calcul d'une matrice de covariance cyclique (U(ak)) 
sur le second moment consiste a choisir deux vecteurs (x(fk.t), x (-f m ,t)) de 

frequence f k et -f m , a les blanchir (5) pour obtenir deux nouveaux vecteurs 
15 (z(^,t),z(-f m ,t)), a calculer K matrices U(t) donnees par la relation 

U(t)=z(f kl t)xz(-f rn ,t) t correspondant a la duree d'observation T, a 
calculer (6) chaque composante u Wak) de la matrice de covariance cyclique 
(U(ak)) en realisant une transformee de Fourier rapide sur K points en 
prenant les K points de meme indice ij des K matrices (U(t)). 

20 

5. Procede selon Tune quelconque des revendications 2 a 4, 
caracterise en ce que pour chaque frequence cyclique (ak) de I'espace 
frequence cyclique - frequence harmonique donne le rapport de 
vraisemblance (V(ak)) est calcule suivant la relation : 

25 V(ok) = -2BT In |det(l - U(<xk). U(ak) + )] 

avec I la matrice identite et B=B canal . 

6. Procede selon I'une quelconque des revendications 1 a 5, 
caracterise en ce que le seuil de detection est determine par une loi du 

30 chL2. 

7. Procede selon I'une quelconque des revendications 1 a 6, 
caracterise en ce que les couples d'antennes (A1.AN) sont tels que les 
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maximums maximorum de sensibility des deux antennes pointent dans deux 
directions orthogonales. 

8. Precede selon Tune quelconque des revendications 1 a 6, 
5 caracterise en ce que les couples d'antennes (A1.AN) sont tels que les 
maximums maximorum de sensibilize des deux antennes pointent dans la 
meme direction. 
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